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Abstract 



For a hypothetical planet on a highly eccentric or- 
bit, we have calculated the osculating orbital pa- 
rameters and its closest approaches to Earth and 
Moon over a period of 750 kyr. The approaches 
which are close enough to influence the climate of 
the Earth form a pattern comparable to that of the 
past climatic changes, as recorded in deep sea sedi- 
ments and polar ice cores. 



1 Motivation 

The information on Earth's past climate obtained 
from deep sea sediments and polar ice cores indi- 
cates that the global temperature oscillated with a 
small amplitude around a constant mean value for 
millions of years until about 3.2 million years ago. 
From then on the mean temperature decreased step- 
wise and the amplitude of the variations increased. 
For the last one million years Earth's climate was 
characterized by a distinct 100 kyr periodicity in 
the advancing and retreating of the polar ice sheets, 
known as glacial-interglacial cycles. This Ice Age 
period ended abruptly 11.5 kyr ago. During the fol- 
lowing time interval, the Holocene, the global mean 
temperature quickly recovered to a value compara- 
ble with that observed 3.2 Myr ago 0, §. Wolfli et 
al. H proposed that during all this time an object 
of planetary size, called Z, existed in a highly ec- 
centric orbit, and showed that Z could be the cause 
of these changes of Earth's climate. Because of the 
postulated small perihelion distance, Z was heated 
up by solar radiation and tidal forces so that it was 
surrounded by a gas cloud. Whenever the Earth 
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crossed this cloud, molecules activating a Green- 
house effect were produced in the upper atmosphere 
in an amount sufficient to transiently enhance the 
mean surface temperature on Earth. Very close 
flyby events even resulted in earthquakes and vol- 
canic activities. In extreme cases, a rotation of the 
entire Earth relative to its rotation axis occured in 
response to the transient strong gravitational inter- 
action. These polar shifts took place with a fre- 
quency of about one in one million years on the 
average. The first of them was responsible for the 
major drop in mean temperature, whereas the last 
one terminated Earth's Pleistocene Ice Age 11.5 kyr 
ago. At present, planet Z does not exist any more. 
The high eccentricity of Z's orbit corresponds to 
a small orbital angular momentum. This could be 
transferred to one of the inner planets during a close 
encounter, so that Z plunged into the sun. Alterna- 
tively, and more likely, Z approached the Earth to 
less than the Roche limit during the last polar shift 
event. In this case it was split into several parts 
which lost material at an accelerated rate because 
of the reduced escape velocity, so that eventually all 
of these fragments evaporated during the Holocene. 

Here, we describe the method used to calculate 
the motion of such a hypothetical object in the pres- 
ence of the other planets of the solar system, and its 
close encounters with the Earth. The calculations 
neglect possible losses of mass, orbital energy and 
angular momentum of Z due to solar irradiation and 
tidal effects. We also disregard the disappearance of 
Z following the last polar shift event. Consequences 
of these effects are discussed in ref. || . 

2 The method 

In order to study the motion of an additional planet 
Z in the gravitational field of the sun and the other 
planets, a set of coupled ordinary differential equa- 
tions (ODE) has to be integrated numerically start- 
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ing at a given time with the known or assumed set 
of orbital parameters of all celestial bodies of in- 
terest. For the calculation we used the Pascal pro- 
gram Odeint which is based on the Bulirsch-Stoer 
method. It includes the modified midpoint algo- 
rithm with adaptive stepsize control and Richard- 
son's deferred approach to the limit 0. The start 
values for the known planets were taken from ref. 
||, where they are given in barycentric rectangu- 
lar coordinates and velocity components referred to 
the mean equator and equinox of J2000.0, the stan- 
dard epoch, corresponding to 2000 January 1.5. All 
values are given in astronomical units (AU) and 
AU/day, respectively. To obtain the heliocentric 
coordinates referred to the ecliptic, we subtracted 
the coordinates of the sun and rotated the resulting 
values for all planets by 23°26'21.448" , the angle 
between the mean equatorial plane and the ecliptic 
of J2000.0. The Earth and the Moon were treated 
separately. The masses of the planets are from ref. 

. For reasons explained in ref. fa] we assume that 
Z was a mars-like object with Mz = 0.11Mb and 
that its hypothetical orbit at the epoch J2000.0 is 
determined by the following heliocentric parame- 
ters: 



semi-major axis: 


a 


= 0.978 


numerical eccentricity: 


e 


= 0.973 


inclination: 


i 


= 0° 


longitude of the perihelion: 


Q 


= 0° 


argument of the perihelion: 




= 0° 


mean anomaly: 


M 


= 270.0° 



The calculation was a classical point-mass inte- 
gration without relativistic corrections. In order 
to save computing time we ignored the influences 
from the three outermost planets Uranus, Neptune 
and Pluto, and restricted the numerical accuracy 
per integration step, the tolerance level eps, to a 
value of 10 -13 Several tests were made to 

check whether these simplifications are acceptable 
or not. First of all, we evaluated the osculating or- 
bital parameters for the Earth over the past 300 kyr 
without Z, but for three different tolerance levels, 
eps = 10~ 13 , 10~ 15 and 10~ 16 , and found that the 
positions of all planets except that of Mercury were 
reproducible to within less than 100 km. A com- 
parison of the Earth's eccentricity and inclination 
variations with the corresponding values published 
by Berger and transferred to the invariable plane 
of the solar system by Quinn et al ^ also showed 
good agreement. We also determined the total an- 



gular momentum of the solar system and found a 
negligible small linear change at the eleventh digit 
of its value. All calculations were performed with 
a 133 MHz PC. Including Z, about 15 h were re- 
quired to cover a time span of 10 kyr with an accu- 
racy of 10 -13 . The sequence of close encounters of 
Z with the inner planets and the Moon amplify the 
errors so that the calculated result represents a pos- 
sible orbit only. Integrations with Z were performed 
both, forward (+300 kyr) and backward (-450 kyr) 
in time relative to J2000.0 in order to demonstrate 
that the behavior of the orbital parameters of Z and 
the Earth were not affected by the change in time 
direction. 

3 Results 

Fig. 1 shows the time dependent variations of the 
osculating orbital parameters of Z (semi-major axis, 
eccentricity, inclination, ascending node and argu- 
ment of the perihelion) over a time period of 750 
kyr. The inclination is defined here as the angle 
between Z's orbital plane and the invariable plane, 
which is perpendicular to the orbital angular mo- 
mentum of the solar system. It is close to the orbital 
plane of Jupiter. The movement of Z is strongly in- 
fluenced by the inner planets and to a lesser extent 
by Jupiter. Of interest is the behavior of the inclina- 
tion, which on average oscillates with a periodicity 
of about 7 kyr only. This period is more than an 
order of magnitude shorter than that of Earth's in- 
clination which is essentially determined by Jupiter 
and, without Z included in the calculation, amounts 
to 100 kyr. The corresponding osculating parame- 
ters for Earth's orbit are displayed in Fig. 2; as men- 
tioned above, they are only marginally disturbed by 
the presence of Z. Sudden jumps in one or several 
orbital parameter values of Z indicate close encoun- 
ters with one of the inner planets. 

In order to find out how often Z approached the 
Earth to distances close enough to influence its cli- 
mate or even to provoke a polar shift ||, we de- 
termined the distance between Z and the Earth for 
each integration step and fitted a parabolic func- 
tion into the values close to the distance of closest 
approach. The minimum of this quadratic function 
was then identified with this distance. The upper 
panel of Fig. 3 shows the result of this evaluation 
for the Z-Earth system over the whole time range 
considered here. Plotted are all encounters with 
distances of less than 0.02 AU = 3 • 10 6 km, as in- 
dicated by the endpoints of each vertical line. The 
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kyr since J2000.0 

Figure 1: Osculating orbital parameters (a, semi-major axis; e, eccentricity; i, inclination; ft, longitude 
of the perihelion; uj, argument of the perihelion) over 750 kyr for the proposed object of planetary size 
Z with mass Mz = 0.11Mb. All parameters are evaluated relative to the invariable plane which is 
perpendicular to the orbital angular momentum of the solar system. The calculation was started at t = 
with the orbital parameter values listed in section 2. Note the "high" frequency of the inclination. 
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Figure 2: Osculating orbital parameters for the Earth (a, semi-major axis; e, eccentricity; i, inclination; 
tt, longitude of the perihelion; u>, argument of the perihelion) for the same time range as in Fig. 1. A 
comparison with the values calculated by Berger 0] and tranformed by Quinn et al. |j] into the invariable 
plane suggests that Earth's orbit is only marginally perturbed by Z. 
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Figure 3: Upper panel: Closest approaches of Z to the Earth over 750 kyr. The plot contains all distances 
of less than 0.02 AU = 3 • 10 6 km, as indicated by the lower endpoints of each vertical line. Lower panel: 
Expanded view over 150 kyr, i.e. from to + 150 kyr, showing details of the irregular clustering of these 
events. This structure sensitively depends on the mass of Z and on the initial orbital parameters selected 
at t = 0. 
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Figure 4: Left side: Closest approaches of Z to the Earth below twice the Moon-Earth distance for 
the same time range. Empty and filled dots correspond to incoming and outgoing movements of Z, 
respectively. Right side: Analogous representation for the approaches to the Moon. The two horizontal 
dashed lines in both figures mark Moon's distance (presently 384 000 km) and the critical distance (30 000 
km), respectively, below which significant polar shifts on Earth as well as dramatic changes in Moon's 
orbital parameters have to be expected. 
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lower panel of Fig. 3 shows details of the irregular 
structure of these encounters which are the result 
of the complex time dependence of the coordinates 
of Z and Earth. Not surprisingly, the encounter fre- 
quency is enhanced whenever the two inclinations 
nearly coincide. Fig. 4 (left side) shows that Z ap- 
proaches the Earth many times to less than Moon's 
distance. As explained in ref. || these flyby events 
can excite strong earthquakes and volcanic activi- 
ties. For distances smaller than 30'000 km such en- 
counters could even result in a rotation of the Earth 
by as much as 20° with respect to the direction 
of the invariant angular momentum. Flyby events 
having distances larger than about the Moon-Earth 
distance are harmless in this respect, but they still 
may influence Earth's climate. In ref. || we have 
shown that Z was surrounded by a gas cloud which 
had an estimated radius of about 2.8 Mio. km at the 
intersection point with Earth's orbit. The interac- 
tion of this cloud with Earth's atmosphere produced 
Greenhouse gases in sufficient amounts to signifi- 
cantly increase the global temperature. 

Close encounters with the Earth also imply close 
encounters with the Moon. These are plotted on the 
right side of Fig. 4. Since the mass of the Moon is 
much smaller than that of Z, recoil effects are much 
larger than in the case of the Earth and, therefore, 
may significantly influence the Moon's orbit. In 
fact, Fig. 5. shows events in which the semi-major 
axis a suddenly changes at some given time by up to 
9% relative to the mean value. Since, according to 
Kepler's third law, the orbital period of the Moon is 
proportional to a 3 / 2 , the orbital angular frequency 
u) ar b also changes, so that the apparent rotational 
frequency Q — uj rot — La or b becomes different from 
zero, assuming that uj ro t, the rotational angular fre- 
quency, is not influenced by such an event. In ref. 
H we propose that the last close encounter of Z 
with the Earth took place only 11.5 kyr ago, so 
that Moon's orbit could also have changed at that 
time. Therefore, the question arises whether the 
tidal friction on the Moon is large enough to stop fl 
within the Holocene. In the appendix we show that 
this is likely to be the case. 

4 Conclusion 

The calculations presented here show that an ob- 
ject of planetary size in a highly eccentric orbit ap- 
proaches the Earth with sufficient frequency to in- 
fluence its climate and even to produce polar shifts, 
the last of which terminated Earth's Ice Age period, 



as explained in ref. 0. Although the pattern of 
these approaches compares well with the observed 
pattern, we have to point out that for various rea- 
sons a one to one correspondence between "theory" 
and "observation" cannot be expected: First, since 
Z no longer exists at present, we have to calculate its 
long-term behavior on the basis of assumed orbital 
parameters and on estimations of its mass. Second, 
Z lost substantial amounts of mass, orbital energy 
and angular momentum each time it passed through 
the perihelion. These effects are neglected in our 
point-mass model calculation. They may signifi- 
cantly alter the orbit of Z, and have to be included 
in an attempt to find out whether the disturbances 
of the orbits of the inner planets due to Z are within 
the boundaries set by present day observations. In 
ref. 3, we also point out that Z disappeared from the 
solar system either during the proposed polar shift 
event 11.5 kyr ago or later on during the Holocene. 
A detailed study of the mechanisms responsible for 
this removal is another important task. 



Appendix: 
Moon's rotation 



Synchronising 



Since a close encounter between Z and the Moon 
can lead to an apparent rotational frequency Q = 
uj ro t — L0 O rb of the Moon which no longer vanishes 
on the average, it is important to know how fast 
tidal friction diminishes f2. The tidal force field is 
parallel to the direction Moon-Earth and has the 
value 

2zGM E R M 



F 



R 3 



(1) 



G is the gravitational constant, Me the mass of the 
Earth, Rm the radius of the Moon, R the distance 
between the centers of Earth and Moon, and z the 
cartesian coordinate in the direction Moon-Earth, 
measured from Moon's center. On Moon's surface 
z = Rm cos(7), where 7 is the angle between the 
z-axis and the direction to the point considered. 
Under the influence of F and the gravitational ac- 
celeration gu on Moon's surface, its shape will de- 
viate from that of a sphere. The deformation is in 
first order given by 



//(-): f/n(i+COs(2 7 )). 



In equilibrium Hq becomes 



GR 2 M M E .. . 
-Wo = - „o = 6.5 m. 



2i?3 5 



M 



(2) 



(3) 
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Figure 5: Moon's orbital changes over 750 kyr induced by the flyby events of Z. Viewed over a sufficiently 
long time (750 kyr) the net effect is small because of the quasi-randomness of these encounters. 



In a time dependent situation elastic tensions re- 
duce the deformation. They will fade away gradu- 
ally so that equilibrium is reached, say, with a re- 
laxation time r. Assuming that the presence of an 
apparent rotation fi results in a deformation which 
lags behind with a phase difference <j> = fir, the 
force field F acting on this deformation will exert 
an angular moment D which tends to stop the rota- 
tion fi. The integration over Moon's surface yields 

sin(20) 16tt G 2 M%R e M p 
2 15" 



D 



R 6 9M 

where p is the density on Moon's surface, 
parent rotation fi varies in time as 



(4) 
The ap- 



dfi 
~dt 



D 



-K 



sin(20) 



(5) 



For the inertial moment of the Moon 2 we use the 
value of a homogeneous sphere, S = ~MmR%, and 
4?r p3 n Then K — m e gr m 



set Mm = 

-16 „-2 



1.0 



3 RmP- - Lllt:11 J1 - M M R 6 

10~ ie s~ 2 . For 4> = fir < 1 the solution of this 
equation is given by 



fi(t) = fi(0) e 



-Krt 



(6) 



The true relaxation time of this deformation is open 
to discussion. If, for example, this time is assumed 
to be r = 1 d = 86 400 s, then the decay constant 
for fi becomes 

rn = = 3500 yr (7) 



Assuming fi = 0.1a;, O i, then the phase becomes <fi = 
0.1oj rot T = 0.023 rad = 1.3°. Thus a rather small 
phase lag of the tidal deformation is sufficient to 
synchronise the lunar rotation within the Holocene. 
This assumed phase lag is smaller than the value 
of 2.16° inferred from astronomical data regarding 
the actual lunar bulge 0. The fact that according 
to Eq. ^ Q(t) never stops, is an artefact of a model 
with a single decay constant. The real dynamics 
also involves slow relaxations, which terminate the 
synchronisation in a finite time. 
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